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METHOD FOR ESTIMATING A SCATTERED RADIATION, PARTICULARLY 
TO CORRECT TOMOGRAPHY OR BOND DENSITOMETRY MEASUREMENTS 

DESCRIPTION 

5 

This invention concerns a method for estimating a 
scattered radiation, for which the main envisaged 
application is the correction of tomography or bone 
densitometry measurements. 

10 The use of an irradiation radiation has the 

disadvantage of producing significant scattered radiation 
through the examined object, especially in the very 
frequent event of a diverging radiation (conical or fan 
shaped) . In other words, each of the detectors situated 

15 behind the object receive not only a primary radiation, 
arising directly from the source by a rectilinear path 
and having crossed through a well defined region of the 
object, but also a scattered radiation from an 
indeterminate source which affects the measurement and 

20 which it would therefore be desirable to correct. 

Several methods are already used. It is in this way 
that the primary radiation can be measured alone if a 
rigorous collimation of the detectors and the source is 
made in order to intercept the scattered radiation, but 

25 in practice said method requires a scanning of the beam, 
which takes time to accomplish, and during which one has 
to put up with movements of the patient, if one is 
examining living organisms. 
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The opposite idea of only measuring the scattered 
radiation has also been proposed. To do this, one 
arranges a discontinuous network of absorbers, such as 
lead balls, between the object and the detectors, in 
5 order to locally stop the primary radiation, in such a 
way that the detectors located behind said absorbers only 
measure the scattered radiation. This process, called the 
"beam stop" process, therefore gives two dimensional 
tables or bundles of the value of scattered radiation, 

10 which one completes by interpolation between the 
detectors located behind the absorbers. The scattered 
radiation estimated in this manner is subtracted from the 
total radiation measured separately. This process is 
precise but has the disadvantage that it imposes two 

15 irradiations of the object and thus the object receives 
double the dose of radiation. A final example of a method 
for correcting the scattered radiation by material means 
comprises the use of anti-scattering grids, but they are 
only partially efficient; it is insufficient for a 

20 conical beam, where the scattered radiation may be 
several times greater than the primary radiation. 

Finally, a certain number of digital methods exist 
for estimating the scattered radiation, from convolutions 
and deconvolutions of measurements, for example: one 

25 could aslo cite French patent 2 759 800 for a different, 
analytical digital method. Said methods are, in general, 
difficult to employ since they depend on parameters 
chosen by the user (convolution kernel, for example) that 
only give good results in favourable situations, such as 



B14164.3/JCI 



3 



small areas where the scattered radiation is low, or 
objects with a relatively homogeneous content. No simple 
method exists that makes it possible, for example, to 
correct the scattered radiation through the thorax or 
5 other major anatomical regions, which are frequently- 
examined but which are unfavourable for correcting the 
scattered radiation due to their very volume and the 
heterogeneity due to the presence of a complex bone 
structure and in which the radiation attenuation capacity 
10 is very different to that of soft tissue. 

Finally, we should mention American patent 6 018 565 
which describes a mixed method, using "beam stop" and 
convolution . 

An essential aim of the invention is to propose a 
15 method for estimating and correcting scattered radiation 
that could be suited to difficult situations of non- 
destructive control of inert or animate objects, or all 
image reconstruction methods. 

The method according to the invention is, in its 
20 most general form, an object imaging method by multiple 
acquisition, comprising an estimation of a scattered 
radiation from an initial radiation that has passed 
through an object while undergoing an attenuation that 
allows a total measurement radiation to pass through, 
25 characterised by: 

- for at least one acquisition, taking a table of 
measurements of a scattered radiation, obtained by 
passing the initial radiation through an imitation of the 
object , 
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and for each of the acquisitions, calculating 
transposition coefficients between the imitation and the 
object, from the initial radiation, the total measurement 
radiation through the object and a total measurement 
5 radiation through the imitation. 

- and weighting of the table of measurements with 
the transposition coefficients. 

Advantageously, the imitation will be a block of 
constant thickness and in a homogeneous material, with an 

10 attenuation similar to a base material of the object; in 
general, the taking of a table of measurements will be a 
selection in a series of tables of measurements of 
scattered radiation, obtained beforehand by successively 
passing the initial radiation through a respective series 

15 of imitations of the object, of different but constant 
thicknesses; and the selection will be made by comparing 
a value of the total measurement radiation through the 
object and a value of the total measurement radiation 
through the imitations. 

20 The weighting coefficients are generally ratios of 

values of a same functional calculated for the object and 
for the imitation. The functional used may be equal to 
the product of the total measurement radiation times the 
logarithm of the ratio of the total measurement radiation 

25 and the initial radiation. 

The invention will now be described in referring to 
the figures, in which: 

- Figure 1 is a general view of an acquisition of 
measurements ; 
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- Figure 2 is a view of a calibration acquisition; 

- and Figure 3 illustrates the steps of the method. 
We will first refer to Figure 1, where an X-ray tube 

1 emits a conical beam 2 towards an object 3 to be 
5 examined (here a patient stretched out on a table 4) 
then, through him, towards a flat network 5 of detectors 
6 arranged in a matrix. The detectors 6 are connected to 
an acquisition device 7 and measure a scattered radiation 
that is superimposed on the primary radiation, only 
10 suitable for the examination or the control of the 
object . 

The estimation of the scattered radiation through 
the patient (object 3) consists firstly in obtaining two 
dimensional tables or bundles of scattered radiation 

15 obtained under comparable circumstances. To achieve this, 
one carries out calibration radiations through imitations 
8 of the object 3 to be examined, in accordance with 
figure 2: the irradiation conditions remain the same, in 
other words one continues to use the tube 1, the beam 2, 

20 the network 5 of detectors 6 and the acquisition device 
7, the imitation 8 however replacing the patient; we have 
also added a grid 9 of lead balls 10 between the 
imitation 8 and the network 5. It ensues from this 
arrangement that the rays 11 going through the balls 10 

25 are completely absorbed and the regions 12 of the network 

5 located in the same line as said rays 11 have detectors 

6 that only measure the scattered radiation at these 
places. One only has to record these measured values and 
interpolate between the regions 12 in order to suitably 
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estimate the scattered radiation from the imitation 8 for 
all of the detectors 6 of the network 5. 

The imitation 8 should be similar to the object in 
order that the radiation scattered by them are identical. 
5 A perfect similarity is not possible, and for this reason 
one makes do with an imitation 8 resembling the object 3 
and for which the associated bundle of scattered 
radiation will be corrected later in order to evaluate 
that of the object. In practice, the imitation 8 may be a 
10 block of a homogeneous material and which has the same 
attenuation coefficient as the base material of the 
object 3: in the case of a human body, essentially 
composed of soft tissue, it is known that Plexiglas 
(polymethacrylate) is suitable. 
15 In order to allow varied measurements, in reality 

one will use several bundles of scattered radiation, 
obtained for as many imitations 8, which will only differ 
by their thickness and therefore by the length of the 
path covered by the rays 11. Said bundles are recorded in 
20 a data base before the useful measurements on the objects 
3, in order to take a bundle of scattered radiation 
comparable to that of an object 3, one will select in 
practice one of the bundles of the data base or, even 
better, a bundle that has been obtained by interpolation 
25 calculations between two of said bundles. The selection 
criterion could be defined by means of a specific ray 13 
ending up in a region 14 of the network 5 and which will 
not pass through either the absorbers 10 of figure 2, or 
through the bone tissues of the patient (or more 
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generally portions of the object 3 in which the 
absorption properties are different to the material of 
the imitation 8) in figure 1. The total radiation, 
primary and scattered, received by the region 14 after 
5 having crossed each imitation 8 will act as index to the 
corresponding scattered radiation table, and the selected 
table will have the index at an identical value to the 
total radiation measured at region 14 through the object 
3 . This all corresponds to going from the state El to the 

10 state E2 in the block diagram of figure 3, that we will 
now begin to explain. 

The continuation of the method essentially consists 
in the correction of the table of scattered radiation 
thus selected to adjust it as best as possible to the 

15 bundle of radiation actually scattered by the object 3. 
To achieve this, one uses all of the information 
available, in other words the total radiation received by 
the detectors 6 beyond the object 3 as the selected 
imitation 8. Said total radiation is designated <Dt, the 

20 scattered radiation Od, the initial radiation from the 
tube Oo and the primary radiation G>, the relation <Dt = <D 
+ Od is met. 

One is then in states E3 and E4 of the block diagram 
of figure 3. Then, one transforms the total radiation 
25 values Ot measured for the object 3 and the selected 
imitation 8 by applying the functionals to them. More 
precisely, it is know in the art that $d is proportional 
to Olog (0/ Oo) ; this relation, which is called the Klein 
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and Nishina law, gives a general appearance of the 
scattered radiation, for lack of its intensity. 

The initial radiation Oo is known; the primary 
radiation 0 is not known f but one agrees to apply this 
5 relation in an approximate manner by replacing it by the 
total radiation Ot, in other words the functional 
employed associates each measured value of the total 
radiation Ot with the calculated value Ot log (Ot / Oo) , 
assumed to be close to the scattered radiation Od at this 
10 location; we have reached states E5 and E6 in the block 
diagram. 

The following step consists in carrying out, for 
each of the detectors 6, the ratio of values given by the 
functional for the object 3 and the imitation 8 selected 

15 according to the formula K= Ot log(Ot/Oo) object/Ot log 
(Ot/Oo) imitation. The weighting coefficients K thus 
obtained will serve to deform the bundle of selected 
scattered radiation in state E2 in order to estimate that 
of the object 3. The results again constitute a two 

20 dimensional table or a matrix having identical dimensions 
to that of the radiation tables since it is associated 
with the network 5 of detectors 6. It is therefore 
possible and advantageous to carry out a spatial digital 
filtering of this matrix by applying a low pass filter 

25 that corrects the coefficients K while only conserving 
the lowest frequencies of their variation and thus making 
them probably closer to reality since the scattered 
radiation varies quite slowly from one point to another. 
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When the table of definitive weighting coefficients, 
designated K' , has been obtained (in state E7) , it serves 
to weight the selected table of scattered radiation 
before state E2, in order to obtain a table of radiation 
5 scattered by the object 3 (state E8, which constitutes 
the sought after estimation) ; the formula applied is <£d 
object = K' <Dd object. These estimated object values Od 
may then be subtracted from the total radiation <Dt 
measured by the detectors 6 to estimate the primary 

10 radiation <D and obtain a more precise radiographic image 
of the obj ect 3 . 

This method applies to simple or multiple energy 
radiations; in the latter case, it is repeated separately 
for each of the energies used. 

15 The functional proposed here is not the only one 

that may be used, and the simpler functional <Dd = k<£ 
(here approaching closer in 3>d = k<Dt) , where k is a 
constant, could thus give good results for estimating <Dd. 

As has been stated, this method is particularly 

20 interesting in tomography or bone densitometry, where an 
in-depth image of the object is reconstructed from a 
multitude of acquisitions (irradiations) taken around the 
object by a mobile network of detectors, then by a 
digital combination of said acquisitions (carried out 

25 using techniques that are not covered by the invention) . 
One of the limitations encountered in practice is the 
excessive dose received by the object; however, we have 
seen that the method according to the invention makes it 
possible to reduce the dose required compared to other 
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methods, especially those that impose a double 
irradiation, to subtract the scattered radiation from the 
total radiation. It is felt that this advantage will be 
marked in linear network detector methods (or with 
superimposed lines in a two dimensional network where one 
obtains one or several cross -sections (two dimensional 
images) through the object. 

When several acquisitions k are undertaken, the 
operating method summarised by figure 3 may be applied to 
each of them, an independent estimation of the scattered 
radiation being made each time. However, situations exist 
where the scattered radiation could be considered 
invariable, particularly for objects with rotational 
symmetry, very frequent in non-destructive control 
methods. Several of the steps of figure 3 then become 
unnecessary: one could thus make do with the selection of 
a single imitation 8 for all of the acquisitions k; if 
necessary, only steps E3 , E5, E7 and E8 involving the 
object need to be repeated for each acquisition. The 
choice of a simplification of the method by the omission 
of certain calibrations could be at the discretion of the 
operator. 
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